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Abstract 

The  low  rate  (C/33)  discharge  data  obtained  from  a  calendar  life  study  were  numerically  analyzed  with  a  single  particle  model.  The  simulation 
showed  that  the  stoichiometric  window  for  the  cathode  shrank  with  capacity  fade.  The  change  of  the  stoichiometric  window  for  the  anode  was 
more  complicated.  The  aged  anode  became  less  charged  when  the  capacity  fade  was  caused  mostly  by  the  loss  of  cyclable  lithium  ions.  The  anode 
would  be  charged  to  a  higher  stoichiometric  number  (or  state  of  charge,  SOC)  when  the  capacity  fade  became  controlled  by  the  loss  of  active 
material  in  the  anode. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Lithium  ion  batteries  have  been  widely  used  in  consumer 
electronics  and  they  are  being  actively  researched  for  automo¬ 
bile,  satellite  and  many  other  applications.  A  few  papers  [1-3] 
are  found  in  the  literature  to  study  numerically  the  capacity  fade 
phenomenon  during  the  life  testing  of  lithium  ion  batteries.  Liaw 
et  al.  [1]  used  a  relatively  simple  equivalent  circuit  model  (ECM) 
to  study  the  capacity  fade  in  a  calendar  life  study.  They  were  able 
to  achieve  high  fidelity  simulation  results  with  a  consistent  set 
of  parameters.  Yet  the  parameters  of  the  ECM  model  mostly  do 
not  correspond  directly  to  the  physical  parameters  of  the  lithium 
ion  batteries,  which  is  a  major  drawback  of  the  method.  Stamps 
et  al.  [2]  studied  the  capacity  fade  of  a  lithium  ion  cell  over  1600 
cycles  using  a  hybrid  estimation  algorithm  combining  elements 
of  both  batch  estimation  and  online  moving  horizon  estimation. 
The  methodology  offers  the  ability  to  monitor  the  long-term  con¬ 
ditions  of  the  battery  while  minimizing  the  noise  attributable  to 
the  cycle-to-cycle  variation.  The  battery  model  was  nevertheless 
oversimplified  in  their  work.  Ramadass  et  al.  [3]  used  a  semi- 
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empirical  model  to  study  the  capacity  fade  of  Sony  18650  cells. 
The  empirical  correlations  for  the  initial  state  of  charge  (SOC) 
and  the  surface  film  resistance  of  the  negative  electrode  were 
established  by  simulating  the  experiment  data  at  selected  cycles. 
These  numerical  studies  were  able  to  provide  useful  information 
about  the  capacity  fade  of  the  cells  in  addition  to  the  experiment 
studies. 

Recently  we  reported  [4]  our  experiment  results  on  a  calen¬ 
dar  life  study.  The  loss  of  active  materials  in  the  electrodes  and 
the  loss  of  cyclable  lithium  ions  were  found  to  be  the  major 
capacity  fade  mechanisms  responsible  for  the  cell  capacity  fade 
at  low  C/33  rate.  However,  the  experiment  provided  little  infor¬ 
mation  about  whether  or  not  the  two  mechanisms  played  an 
equally  important  role  throughout  cell  capacity  fade.  In  the  work 
reported  here,  the  experimental  data  were  numerically  analyzed 
to  answer  the  question.  Furthermore,  the  change  of  stoichiomet¬ 
ric  windows  of  the  individual  electrodes  with  capacity  fade  is 
discussed  in  details  based  on  our  simulation  results. 

2.  Experimental 

The  experiment  details  of  the  calendar  life  study  were 
revealed  in  a  previous  paper  [4].  The  test  procedures  were  sum¬ 
marized  in  Fig.  1 .  The  procedures  included  rate  capability  tests, 
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Fig.  1.  Schematic  diagram  of  the  test  procedures  for  the  calendar  life  study  (C  =  1.828  A).  The  cells  were  float  charged  during  storage. 


Nomenclature 

A  electro-active  surface  area  (cm2) 

Cq  Li+  concentration  in  the  electrolyte  (mol  cm-3) 
Cmax  maximum  Li+  concentration  in  the  particles 
(mol  cm-3) 

cs  Li+  concentration  in  the  solid  particles 
(mol  cm-3) 

cs  volume  averaged  Li+  concentration  inside  spher¬ 
ical  particles  (mol  cm-3) 

cse  Li+  concentration  at  the  surface  of  the  solid  par¬ 
ticles  (mol  cm-3) 

Ds  solid  phase  Li+  diffusion  coefficient  in  the  parti¬ 
cles  (cm2  s-1) 

4pp  applied  current  (A) 

7int  intercalation  current  density  (A  cm-2) 
k  kinetic  rate  constant  (A  cm_2)/(mol  cm-3)1,5 

n  number  of  terms  in  Eq.  (5) 

Rs  radius  of  the  spherical  particles  (cm) 

U  open  circuit  potential  of  the  electrode  (V) 

V  volume  of  the  electrode  (cm3) 

vo,/  the  stoichiometric  number  of  the  electrode  at 
charged  state 

Xi  the  stoichiometric  number  of  the  electrode  at  dis¬ 
charged  state 

Greek  letters 

ss  volume  fraction  of  the  active  material  in  the  elec¬ 
trode 

Xm  positive  eigen-value  determined  from 
Xm  —  tan  (A,m). 

0  electrode  potential  (V) 


storage  at  float  charge  condition  (4.1V)  and  capacity  checks. 
The  experiments  were  conducted  at  five  temperatures:  5, 15,  25, 
35  and  45  °C.  The  capacity  fade  measured  at  the  C/33  rate  is  pre¬ 
sented  in  Fig.  2.  The  capacity  fade  at  5  and  15  °C  is  linear  with 
time,  while  the  nonlinear  capacity  fade  at  25,  35  and  45  °C  might 
also  be  viewed  as  two  linear  capacity  fade  regimes  in  sequence. 


Our  half-cell  study  [4]  showed  that  the  used  electrodes  har¬ 
vested  from  the  aged  cells  had  reduced  capacities  (discharged  at 
very  low  rate)  as  compared  to  fresh  electrodes,  indicating  a  loss 
of  active  materials  in  the  electrodes  with  cell  capacity  fade.  Our 
conclusion  [4]  is  that  the  loss  of  active  materials  in  the  electrodes 
and  the  loss  of  cyclable  lithium  ions  are  the  two  major  capacity 
fade  mechanisms  responsible  for  the  cell  capacity  fade  at  low 
C/33  rate. 

Nevertheless,  the  experiment  does  not  supplied  answers  to 
the  question  that  whether  or  not  the  two  mechanisms  played 
an  equally  important  role  throughout  cell  capacity  fade,  such 
as  the  nonlinear  capacity  fade  exhibited  at  high  temperatures. 
Actually  it  is  quite  difficult  to  answer  the  question  experimen¬ 
tally,  because  a  large  sample  size  of  pouch  cells  and  regular 
destructive  analysis  are  required  for  that  purpose.  That  is,  we 
need  to  open  a  few  pouch  cells  (thus  useless  for  further  aging 
test)  regularly  during  the  aging  study  in  order  to  experimen¬ 
tally  measure  the  severeness  of  active  material  loss  in  individual 
electrodes. 

In  this  work,  an  alternative  numerical  approach  was  applied 
to  investigate  the  roles  of  the  two  capacity  fade  mechanisms 


Fig.  2.  The  capacity  fade  of  the  cells  in  the  calendar  life  study. 
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at  different  temperatures.  In  the  following  section,  a  so-called 
single  particle  model  [5,6]  was  introduced,  and  then  used  to 
analyze  the  low  rate  (C/33)  discharge  data.  The  results  from  our 
simulation  were  discussed  in  details  in  Section  4. 


3.  Model 


The  single  particle  model  [5,6]  uses  the  assumption  that  the 
behavior  of  the  porous  electrodes  can  be  represented  by  spherical 
particles  of  the  active  materials  in  the  electrodes.  The  concentra¬ 
tion  and  potential  distributions  in  the  solution  phase  are  assumed 
to  be  negligible,  which  is  a  reasonable  assumption  under  low  rate 
situations.  The  model  equations  for  the  single  particle  model 
consist  of  the  solid  phase  diffusion  equation  (Fick’s  diffusion 
law)  in  a  spherical  particle  and  the  Butler  Volmer  equation  for 
each  electrode: 


9cs  i  19/  9  9cs  i 

-fr  =  Z2^(D^2^r 


dt  r 2  dr 
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The  intercalation  current  density  j[ntj  is  calculated  based  on 
the  equivalent  electro- active  surface  area  (A;)  of  the  electrode 


■  _  ^app 

7int,z  —  . 

Ai 


((3/J?M)fis.,-Vi-) 


(3) 


where  V*  is  the  volume  of  the  electrode  calculated  based  on  the 
length,  width  and  thickness  of  the  electrode.  It  is  assumed  to 
be  constant  during  the  aging  process.  ssj  is  the  volume  fraction 
of  the  active  material,  which  is  defined  as  the  ratio  between  the 
volume  of  the  active  material  over  the  volume  of  the  electrode 
Vi.  The  decrease  of  the  volume  fraction  ssj  indicates  a  loss  of 
active  material  in  the  electrode,  which  could  be  caused  by  the 
isolation  of  active  material  particles  from  the  electrode  matrix 
(mechanical  failure),  etc. 

The  pseudo  steady  state  method  (PSS)  [7,8]  was  used  to  sim¬ 
plify  the  diffusion  equation  from  a  partial  differential  equation 
(PDE)  with  time  and  space  dependence  to  a  set  of  differential 
and  algebraic  equations  (DAEs)  with  only  time  dependence. 
The  PSS  equations  for  the  diffusion  in  the  solid  particles 
are 


_  3  Jint,z 

dt  Rsj  F 


(  I  _  7int,z  Rs,i  2  Rs  i  Jint,z 
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where  the  number  of  terms  4V’  in  the  summation  is  typically 
set  equal  to  5;  csj  is  the  volume  averaged  concentration  in  the 
spherical  particle  and  cs?z  |r=/?  is  the  concentration  at  the  sur¬ 
face  of  the  particle.  The  state  of  charge  or  the  stoichiometric 
number  of  the  electrode  is  defined  as 


£max,z 

The  km  values  are  determined  from  the  eigen- value  equation: 


7-m  —  t^n  (Xm ) 


(7) 


The  single  particle  model  is  used  to  fit  the  low  rate  (C/33) 
experimental  data  from  the  rate  capability  tests.  Four  model 
parameters  are  found  to  change  with  storage  time  when  fitting 
the  C/33  rate  discharge  curves.  The  model  parameters  are  initial 
state  of  charge  (vo)  of  the  anode  and  the  cathode  at  the  beginning 
of  discharge  and  the  volume  fraction  (es)  of  the  active  materials 
in  both  electrodes.  The  least  square  regression  method  [9]  was 
used  to  fit  the  experimental  data  with  the  single  particle  model. 
It  was  assumed  that  the  kinetic  and  transport  limitations  are  not 
important  during  the  low  rate  discharge  in  the  fitting  process. 


4.  Results  and  discussion 


Fig.  3  compares  the  C/33  discharge  profiles  for  fresh  cells 
stored  at  different  temperatures.  The  parameter  values  are  com¬ 
pared  in  Fig.  4.  The  values  of  the  volume  fractions  (es)  are  close, 
indicating  that  the  fresh  cells  have  similar  amounts  of  active 
materials  in  the  electrodes.  The  SOC  of  the  LiCoC>2  electrode 


Fig.  3.  Discharge  profiles  for  fresh  cells  stored  at  different  temperatures 
(C=  1 .656  A).  The  profiles  from  the  top  to  the  bottom  are  45-5  °C.  The  points 
are  experiment  data  and  the  lines  are  the  simulated  profiles. 
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Fig.  4.  The  fitting  parameter  values  for  fresh  cells  at  different  temperatures. 
Both  the  LiCoC>2  and  carbon  electrodes  become  less  charged  (i.e.,  larger  xo,p0s 
and  smaller  xo,neg)  at  low  temperatures. 


Capacity  (Ah) 


Fig.  6.  The  experimental  discharge  profiles  taken  in  the  rate  capability  tests 
at  25  °C.  The  fact  that  less  ripples  are  observed  at  the  beginning  of  discharge 
profiles  indicates  that  the  LiCoC>2  electrode  is  gradually  less  charged. 


(*0,pos)  at  the  beginning  of  discharge  decreases  as  the  temper¬ 
ature  increases.  The  cathode  becomes  more  charged  at  higher 
temperatures.  This  is  directly  reflected  in  the  experimental  pro¬ 
files  as  there  are  fewer  ripples  on  the  top  of  the  discharge  profiles 
at  low  temperatures,  as  shown  in  Fig.  3.  The  ripples  are  char¬ 
acteristic  of  the  open  circuit  potential  (OCP)  of  the  LiCoC>2 
electrode,  as  shown  in  Fig.  5.  The  presence  of  these  characteris¬ 
tic  ripples  makes  it  relatively  easy  to  locate  vo,pos  for  the  LiCo02 
electrode.  By  carefully  observing  the  ripples  on  the  low  rate  pro¬ 
files,  one  can  tell  whether  the  LiCoC>2  electrode  is  charged  to  a 
higher  or  lower  stoichiometric  number  (or  SOC).  For  example, 
Fig.  6  shows  the  C/33  rate  discharge  profiles  measured  at  25  °C 
in  the  rate  capability  tests.  The  experimental  data  show  that  there 
are  fewer  ripples  at  the  beginning  of  the  discharge  profiles  after 


storage.  This  indicates  that  the  aged  LiCoC>2  electrodes  become 
increasingly  less  charged  after  storage,  even  though  the  end  of 
charge  voltage  (EOCV)  for  the  full  cell  remained  the  same.  The 
simulation  also  shows  that  the  carbon  electrode  becomes  less 
charged  at  low  temperatures  (with  smaller  *o,neg)>  which  causes 
the  cells  to  deliver  less  capacity. 

Figs.  7  and  8  present  comparison  of  the  experimental  and 
simulated  discharge  profiles  at  C/33  rate  for  5  and  15  °C,  respec¬ 
tively.  The  changes  of  the  four  fitting  parameters  with  storage 
time  are  displayed  in  Figs.  9  and  10,  respectively.  The  capacity 
fade  and  the  change  of  fitting  parameters  at  5  and  15  °C  follows 
similar  pattern.  The  aged  LiCoC>2  cathode  becomes  increas¬ 
ingly  less  charged,  as  the  stoichiometric  number  (vo,pos)  °f  the 
cathode  increases  slightly  with  time.  The  anode  also  becomes 
less  charged  (with  smaller  vo,neg)«  The  volume  fraction  of  active 
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Fig.  5.  The  open  circuit  potentials  for  the  LiCoC>2  and  carbon  electrode. 


Fig.  7.  Comparison  of  experiment  and  simulated  discharge  profiles  at  5  °C. 


Parameter  Value  Parameter  Value  Voltage  (V) 
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8.  Comparison  of  experiment  and  simulated  discharge  profiles  at  15  °C. 


Fig.  9.  Change  of  fitting  parameter  values  with  storage  time  at  5  °C. 


Fig.  10.  Change  of  fitting  parameter  values  with  storage  time  at  15  °C. 


Fig.  11.  Comparison  of  experiment  and  simulated  discharge  profiles  at  35  °C. 

material  in  the  anode  £s,neg  decreases  after  storage  at  float  charge 
condition.  This  indicates  that  some  electro-active  surface  area 
becomes  inaccessible  after  storage  because  of  the  deposition 
of  side  reaction  products  and/or  some  carbon  particles  become 
mechanically  isolated  from  the  electrode  and  inactive  for  the 
insertion  reaction.  The  parameter  £s?pos  does  not  change  much 
with  storage  at  these  temperatures.  The  model  parameters  are 
shown  to  change  largely  in  a  linear  pattern,  consistent  with  the 
linear  capacity  fade  of  the  full  cells. 

Figs.  1 1  and  12  compare  the  experimental  and  simulated  dis¬ 
charge  profiles  at  35  and  45  °C.  The  changes  of  fitting  parameter 
values  with  storage  time  are  shown  in  Figs.  13  and  14.  The 
volume  fractions  drop  slightly  in  the  first  few  months  of  stor¬ 
age  where  the  capacity  fade  is  small.  Then  the  volume  fraction 
drops  at  a  much  faster  pace  and  the  capacity  fade  of  the  cells 
also  accelerates.  The  decrease  of  the  volume  fractions  of  active 


Fig.  12.  Comparison  of  experiment  and  simulated  discharge  profiles  at  45  °C. 
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Fig.  13.  Change  of  fitting  parameter  values  with  storage  time  at  35  °C. 


materials  is  much  more  significant  at  high  temperatures  than  that 
at  low  temperatures.  These  results  suggest  that  there  could  be 
two  different  mechanisms  controlling  the  capacity  fade  at  these 
temperatures.  The  capacity  fade  was  initially  caused  mostly  by 
the  loss  of  lithium  ions  due  to  side  reactions.  Then  the  capacity 
fade  becomes  controlled  by  the  rapid  loss  of  active  material  in 
the  electrodes.  Our  half-cell  study  [4]  showed  that  there  were 
significant  losses  of  active  material  in  the  used  electrodes  taken 
from  the  aged  cells.  We  think  the  loss  is  mostly  caused  by  the 
mechanical  changes  of  the  electrodes,  where  the  solid  particles 
gradually  lose  contact  with  the  electrode  during  storage  (at  float 
charge  condition).  The  stoichiometric  number  for  the  anode  in 
the  charged  state  (vo,neg)  decreases  when  the  capacity  fade  is 
mostly  caused  by  the  loss  of  lithium  inventory.  However,  when 
the  capacity  fade  is  controlled  by  the  loss  of  active  materials, 
the  vo,neg  starts  to  increase,  indicating  that  the  anode  is  charged 


Fig.  14.  Change  of  fitting  parameter  values  with  storage  time  at  45  °C. 


Fig.  15.  Comparison  of  experiment  and  simulated  discharge  profiles  at  25  °C. 

to  a  higher  SOC.  A  plausible  explanation  for  the  behavior  of 
■*o,neg  is  that  the  active  material  loss  in  the  anode  is  so  severe 
that  it  outpaces  the  loss  of  lithium  inventory  in  the  system. 
Thus,  the  remaining  active  carbon  material  in  the  anode  needs 
to  be  charged  to  a  higher  SOC  to  accommodate  the  lithium  ions 
transferred  from  the  cathode.  The  stoichiometric  number  of  the 
cathode  gradually  increases  during  capacity  fade,  meaning  the 
cathode  becomes  less  charged.  This  result  is  consistent  with 
the  experimental  observation;  that  is,  the  characteristic  ripples 
from  the  cathode  gradually  disappear  in  the  discharge  profiles 
as  storage  time  increases. 

Comparisons  of  the  experimental  and  simulated  profiles  at 
25  °C  are  presented  in  Fig.  15.  The  changes  of  fitting  parameters 
with  storage  time  are  shown  in  Fig.  16.  The  model  parame¬ 
ters  change  in  similar  patterns  as  those  at  higher  temperatures, 
except  the  stoichiometric  number  of  the  anode  (vo,neg)-  The  vo,neg 


Fig.  16.  Change  of  fitting  parameter  values  with  storage  time  at  25  °C. 
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Fig.  17.  The  simulatedxpos  (SOC  at  the  end  of  low  rate  discharge)  for  the  LiCoC>2 
electrodes  after  storage. 

decreases  continuously  with  capacity  fade  at  25  °C.  An  expla¬ 
nation  for  this  could  be  that  the  loss  of  active  carbon  material  at 
25  °C  is  not  severe  enough  to  outpace  the  loss  of  lithium  inven¬ 
tory.  Thus  loss  of  active  material  in  the  anode  could  not  trigger 
the  change  of  the  parameter  vo,neg  at  25  °C  as  it  does  at  higher 
temperatures. 

With  the  model  parameters  readily  available  for  each  dis¬ 
charge  curve,  it  is  of  interest  to  see  what  the  stoichiometric 
numbers  (SOC)  would  be  at  the  end  of  the  low  rate  discharge 
for  the  cathode  and  the  anode. 

Fig.  17  shows  the  simulated  stoichiometric  numbers  for  the 
cathodes  stored  at  different  temperatures.  When  cyclable  lithium 
ions  are  irreversibly  lost  in  the  charge/discharge  process,  there 
are  not  enough  of  them  available  in  the  cell  to  intercalate  the 
cathode  fully  back  to  the  fresh  state.  The  consequence  of  this 
is  that  the  cathode  is  gradually  emptied,  as  confirmed  by  the 
increasingly  smaller  stoichiometric  numbers.  The  smaller  the 
stoichiometric  number  (vpos  in  the  figure),  the  less  the  LiCoC>2 
electrode  is  intercalated.  Fig.  17  also  shows  the  temperature 
effect  on  the  xvos.  As  the  cells  lose  capacity  (lithium  inventory) 
faster  and  greater  at  increased  temperatures,  the  stoichiometric 
numbers  also  decrease  faster  and  steeper.  Fig.  18  shows  the  sto¬ 
ichiometric  numbers  for  the  carbon  electrode  (vneg)  at  the  end 
of  discharge  after  storage.  The  carbon  electrode  appears  to  be 
mostly  discharged  to  the  same  vneg  all  the  time  except  at  high 
temperatures,  where  the  anode  becomes  slightly  less  discharged. 

Based  on  our  simulations,  it  is  clearly  shown  that  the  stoichio¬ 
metric  window  of  the  LiCo02  electrode  shrinks  with  capacity 
fade  after  cells  are  stored  for  some  time  under  float  charge  con¬ 
ditions.  That  is,  the  fresh  LiCo02  electrode  might  be  cycled 
between  the  x  marker  (at  charged  state)  and  the  O  marker  (at 
discharged  state)  as  shown  in  Fig.  19  when  the  cell  is  cycled 
between  3.3  and  4.1V.  After  storage,  the  aged  LiCoC>2  elec¬ 
trode  could  not  be  charged  to  the  same  stoichiometric  number 
(x  marker)  even  though  the  charge  protocol  remains  the  same, 
because  of  the  increased  cell  impedance.  The  aged  UC0O2  elec¬ 
trode  could  not  be  discharged  to  the  same  SOC  (Q  marker)  with 


Fig.  18.  The  simulated  xneg  (SOC  at  the  end  of  low  rate  discharge)  for  the  carbon 
electrodes  after  storage. 

the  same  discharge  protocol  because  of  the  loss  of  lithium  inven¬ 
tory.  Consequently,  the  stoichiometric  window  on  the  LiCo02 
electrode  shrinks  with  capacity  fade  in  this  calendar  life  study. 
The  simulation  also  shows  that  the  downward  movement  of  the 
x  marker  is  relatively  small,  but  the  upward  movement  of  the 
O  marker  is  more  obvious  and  is  closely  related  to  the  degree 
of  cell  capacity  fade. 

The  change  of  the  stoichiometric  window  on  the  carbon  elec¬ 
trode  is  more  complicated.  The  fresh  carbon  electrode  might  be 
cycled  between  the  x  marker  (at  charged  state)  and  the  O  marker 
(at  discharged  state)  as  shown  in  Fig.  20.  The  movement  of  the 
x  marker  for  an  aged  carbon  electrode  depends  on  the  control¬ 
ling  capacity  fade  mechanisms  (loss  of  Li+  ions  or  loss  of  active 
carbon).  At  low  temperatures,  the  capacity  fade  is  mostly  caused 
by  loss  of  cyclable  lithium  ions.  The  carbon  electrode  gradually 
becomes  less  charged,  that  is,  the  x  marker  moves  to  the  left 
in  Fig.  20.  At  high  temperatures  (in  Figs.  13  and  14),  when  the 


Fig.  19.  The  shrinkage  of  the  stoichiometric  window  for  the  cathode  with  capac¬ 
ity  fade. 
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Fig.  20.  The  change  of  the  stoichiometric  window  for  the  anode  after  storage 
(at  float  charge  condition).  The  anode  becomes  less  charged  when  the  capacity 
fade  is  caused  mostly  by  the  loss  of  cyclable  lithium  ions.  The  anode  would  be 
charged  to  a  higher  SOC  when  the  loss  of  active  materials  controls  the  capacity 
fade  and  outpace  the  loss  of  lithium  ions. 

capacity  fade  is  initially  controlled  by  the  loss  of  cyclable  lithium 
ions,  the  x  marker  moves  to  the  left  (the  anode  is  less  charged). 
When  the  loss  of  active  materials  becomes  severe  enough  to 
outpace  the  loss  of  lithium  inventory  and  controls  the  capacity 
fade,  the  anode  turns  to  be  charged  to  a  higher  SOC  and  the  x 
marker  moves  to  the  fight  in  Fig.  20.  The  aged  carbon  electrodes 
appear  to  be  discharged  mostly  to  the  same  SOC  (the  O  marker). 
At  high  temperatures,  the  anode  tends  to  become  slightly  less 
discharged  and  the  Q  marker  moves  a  little  bit  downward  in 
Fig.  20. 

5.  Conclusions 

Our  experiment  results  showed  that  the  loss  of  active  mate¬ 
rials  in  the  electrodes  and  the  loss  of  cyclable  lithium  ions  were 
the  major  capacity  fade  mechanisms  responsible  for  the  cell 
capacity  fade  at  low  C/33  rate  in  a  calendar  life  study.  But  the 
experiment  provided  little  information  on  how  the  two  mecha¬ 
nisms  evolve  with  cell  capacity  fade.  In  this  work,  we  analyzed 


the  experimental  data  numerically  with  a  single  particle  model. 
Our  simulation  results  suggest  that  the  cell  capacity  fade  could 
have  different  controlling  mechanisms  at  different  temperatures. 
The  capacity  fade  at  low  temperature  is  mostly  caused  by  the 
loss  of  the  lithium  ion  inventory,  with  only  minor  loss  of  active 
materials  in  the  electrodes.  At  high  temperatures,  the  cell  capac¬ 
ity  fade  exhibits  a  two  regime  pattern,  where  the  loss  of  lithium 
ion  inventory  and  the  severe  loss  of  active  materials  controls  the 
capacity  fade  at  the  respective  regime. 

In  addition,  we  also  investigated  the  changes  of  the  stoichio¬ 
metric  windows  for  the  cathode  and  the  anode  with  cell  capacity 
fade.  The  stoichiometric  window  for  the  LiCoC>2  electrode 
shrinks  with  cell  capacity  fade.  The  change  of  the  stoichiomet¬ 
ric  window  on  the  anode  is  more  complicated  and  depends  on 
the  controlling  capacity  fade  mechanisms.  When  the  capacity 
fade  is  caused  mostly  by  Li+  loss,  the  anode  becomes  consis¬ 
tently  less  charged.  When  the  loss  of  active  material  in  the  anode 
controls  the  capacity  fade  and  outpaces  the  loss  of  lithium  ion 
inventory,  the  anode  is  gradually  charged  to  a  higher  SOC,  in 
danger  of  becoming  overcharged  at  high  temperatures  such  as 
35  and  45  °C. 
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